Shaking the entropy out of a lattice: atomic filtering by vibrational excitations 
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We present a simple and efficient scheme to reduce atom-number fluctuations in optical lattices. 
The interaction-energy difference for atoms in difi^erent vibrational states is used to remove excess 
atomic occupation. The remaining vacant sites are then filled with atoms by merging adjacent wells, 
for which we implement a protocol that circumvents the constraints of unitarity. The preparation of 
large regions with precisely one atom per lattice site is discussed for both bosons and fermions. The 
resulting low-entropy Mott-insulating states may serve as high-fidelity register states for quantum 
computing and as a starting point for investigations of many-body physics. 
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I. INTRODUCTION 

An optical lattice with single-atom occupancy of each 
site is the ideal low-entropy starting point for quantum 
computation [iMa and for the simulation of condensed 
matter systems [a 0] ■ Even at the lowest temperatures 
achieved to date, however, defects, i.e. vacant or mul- 
tiply occupied sites, are still common [8-10] and jeop- 
ardize the aforementioned applications: For these, ex- 
ceptionally low temperatures are required, and neither 
vacancies nor multiple occupation of lattice sites can be 
tolerated. Cooling and thermometry schemes beyond the 
nanokelvin range are therefore being developed [ll[ . The 
purification of the state is a challenge: For an ideal low- 
entropy Mott insulator with single-atom occupancy of 
lattice sites, one needs to reliably eliminate multiple oc- 
cupations using number filtering, and, subsequently, fill 
vacant lattice sites with atoms from a reservoir. 

Multiple occupations - or, for fermions, occupation of 
higher excited states - can be dealt with by global 
manipulation [l3l - [l5| . using the difference in the inter- 
action energy between internal states to make transition 
frequencies density-dependent. In this way higher single 
well occupation numbers can be filtered to unity. 

To avoid vacant sites, one can use single-site resolution 
to first measure the occupation of the lattice in a site- 
resolved manner and to subsequently rearrange atoms 
using spin-dependent lattices, site-selective spin-flips and 
spin-selective shifts [l^. A crucial requirement of this 
proposal is the ability to recool atoms to their vibra- 
tional ground state, which despite recent advances [13] 
still remains an unsolved challenge. Alternatively, one 
can utilize global many-body dynamics [l^ |3 13 where 
entropy is removed followed by one or more periods of 
global equilibration. 

While the combination of filtering and tunneling op- 
erations ideally prepares a lattice with unit occupancy 
of each well, the difference in interaction energy between 
different hypcrfine states may be too small [20[ and the 
thermalization may be too slow for the procedure to com- 
pete with the heating under realistic experimental con- 



ditions. On the other hand, the dependence of the in- 
teraction energy between atoms on their spatial overlap 
implies a considerable difference in the interaction en- 
ergy for different vibrational states. Experimentally, the 
removal of excess atoms was recently performed using 
parametric lattice modulation by exploiting different vi- 
brational states [l3]- The problem of vacancies has not 
been tackled experimentally so far. 

In this Article we introduce a combined purification 
protocol for number-filtering and filling of vacant sites, 
for a thermal Mott insulator. We achieve occupation- 
number-dependent transfer that does not rely on the usu- 
ally small natural difference in interaction strength for 
different hyperfine states, but instead we use the interac- 
tion shifts of different vibrational states [loj . 

We will use these interaction shifts for (i) number fil- 
tering, where we remove excess atoms from the wells (Sec- 
tion |TT| and for (ii) filling of vacancies by merging, where 
we merge each lattice well with two auxiliary traps in a 
unitary way to establish unit filling of each well (Sec- 
tion imp . We then apply our results to the experimen- 
tally relevant case of a shallow external confinement with 
a non-uniform thermal distribution of atoms and vacan- 
cies, in Section HVl We analyze how the central domain 
approaches the perfect state, and how the size of this per- 
fect domain grows with successive merging steps. Section 
fVl summarizes our results. 



II. REMOVING EXCESS ATOMS: 
FILTERING 
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As starting point we consider a system where only a 
particular hyperfine state (a) and the vibrational ground 
state are populated, but with fluctuating occupation 
from site to site. In our numerical examples below, the 
starting point will be low temperature approximations to 
the Mott-Insulating (MI) state where the vacancy prob- 
ability is a function of the interaction strength, the tem- 
perature and the chemical potential (see Eq. ^). 

To filter the population on all sites so that any higher 
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state = 2, |2)g, has a different interaction energy 



FIG. 1: (color online) Filtering step, where the quantum state 
with n = 3 vibrational ground state atoms in (a) is reso- 
nantly coupled by Raman laser beams to the state with one 
vibrationally excited atom in (/3). This atom is subsequently 
removed. The sequence is then repeated for n — 2. 



occupancy is reduced to unity, we propose to drive a Ra- 
man transition from the vibrational ground state of one 
hyperfine state (a) to the second vibrationally excited 
level of another hyperfine state (/3). Since the atomic 
population is transferred between states with different 
spatial distributions (see Fig. [I}, the interaction energy 
in the initial and final state depends significantly on the 
occupation of the states. A pair of atoms residing in 
vibrational states {v) and have an interaction energy 



U^,ii = (2 - S^^fj,)- 
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where the constant of proportionality is related to the s- 
wave scattering length Ogc and includes a factor account- 
ing for exchange |2l| . For simplicity, we ignore in the 
following the difference between the internal state inter- 
action strengths (for ^''Rb C/q/j ~ 0.98Uaa [l^l), and we 
assume anisotropic harmonic oscillator wave functions for 
the lowest vibrational states. In the harmonic approxi- 
mation we find the following relative coupling strengths 
for the lowest levels (states are labeled by the ID vibra- 
tional quantum number in the direction of excitation), 
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Since these interaction strengths differ significantly more 
than the ones of different internal states, much higher 
sensitivity to the distribution of atoms in different vibra- 
tional states is achieved, which allows a faster population- 
dependent transfer. More fundamentally, however, the 
difference enables the algorithmic manipulation discussed 
in the next section. 

The configuration with n atoms in the vibrational 
ground state ~ and hyperfine state (a), |0)^, has 
the interaction energy 
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while the state with n — 1 atoms in |0)^, and 1 atom in 
hyperfine state (/?) and the second vibrationally excited 
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Raman laser beams with the corresponding two-photon 
detuning and a suitable wave number difference can 
transfer exactly one atom to |2)^ in all sites with pre- 
cisely n atoms, while leaving the population unchanged 
in all other sites. This transfer step is illustrated in Fig.[l] 
for rt = 3. Sweeping a sequence of pulses resonant with 
the transfer from traps with rimax, n^ax — 1, • ■ • , 2 atoms, 
each followed by the removal of the atoms in the hyper- 
fine state (/3), will transform sites with 1 < n < rimax 
into sites with exactly one atom in |0)^. Choosing n„iax 
large enough, we are left with defects only in the initially 
vacant sites. 



III. FILLING VACANT SITES BY MERGING 
A. Constraints of unitarity 

To fix the remaining vacancies, we propose to merge 
each well of the filtered system with auxiliary wells and 
to perform operations on the combined wells that result 
in final states with the unit occupancy component sig- 
nificantly increased. We envision this being performed 
in parallel across the entire system by merging planes or 
strings of sites. 

Intuitively, an "OR" -operation on adjacent wells max- 
imizes the probability to find at least one atom in the 
ground state of the desired target well. Given an ini- 
tial probability for vacancies e, merging k wells into one 
could yield a final vacancy probability e*^ , i.e. a dramatic 
improvement can be achieved. Such an "OR" -operation, 
however, is irreversible, and cannot be achieved by uni- 
tary dynamics. For example, in order to merge two wells 
to a target well (here: the left well), the desired mapping 
should read: 
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(5) 



which is clearly non-unitary. At first sight, it is thus not 
evident how an improvement of the vacancy probability 
can be achieved without recourse to read-out of the well 
populations, as proposed in Ref. [l6j . 

Here, we present a unitary protocol that is inspired by 
the principal idea of the "OR" -operation, while it circum- 
vents the apparent constraint of unitarity by exploiting 
the larger state-space of three wells and by using particles 
in one of the wells as ancilla. 

We assume that the vacancy probability in each well 
is e, and show the eight different possible occupations of 
three wells in Tab. HI The middle well is taken to be the 
target well, which is occupied (or not) before the proto- 
col according to nj^. The occupation of the target well 



|*initial) P 


n\n 


nvt 


f 

Mm 




1 


1 


1 


|1,1,0) e(l-e)2 


1 


1 


1 


|0,1,1) 


1 


1 


1 


|1,0,1) 





1 


1 


|0,0,1) e'(l-e) 





1 





|0,1,0) 


1 


1 


1 


|1,0,0) 





1 





|0,0,0) 












(I) 



TABLE I: Overview over probabilities for vacancies in the 
three-well system. The initial occupations of the three wells 
correspond to the Fock-state \ni,nm, rir), the respective prob- 
abilities are then given by P. When the target well is (not) 
occupied initially, n\^-^ = 1 (0). After the application of a 
hypothetical non- unitary "OR" -operation, the target well oc- 
cupation is n'^ , while our protocol gives n^^. 



changes after the hypothetical OR- merger (n^^ ) and af- 
ter our protocol (n^J. As we will show in the follow- 
ing, the final occupation can indeed be achieved via 
unitary operations. Our protocol improves only the oc- 
cupation of the middle well for the initial state |1,0, 1), 
i.e. it requires both neighboring wells to be occupied for 
success. Given the population of this initial state of the 
order of e, the procedure leads to a significant decrease 
of the probability of vacancies. 



B. Merging protocol 

1. Merging of the central and right well 

In a first step, we merge the central target well with 
its right neighboring well, such that the middle (right) 
well ground state is mapped to the lowest (first excited) 
level of the combined well. The process is depicted in 
Fig. [21 where the third (left) well is omitted for clarity, 
since this first part of the protocol leaves it untouched. 

The merging procedure has to preserve the single- 
particle states while at the same time being insensitive 
to the possible interaction between the atoms. Minimiz- 
ing the error in the merging process is complicated but in 
[22| we have numerically optimized the merging efficiency 
and for appropriate parameters show that error probabil- 
ities e„i ~ 10"'' can be obtained in less than 100/is even 
in the presence of interactions. The configuration with 
two atoms (upper line) can be dealt with by transfer- 
ring the excited vibrational state population into another 
hyperfine and vibrational state using an upper sideband 
{5v = -f 1) Raman transition. In principle, this operation 
would also affect the singly occupied instances (lines (II) 
and (III) in Fig. (2); however, since the doubly occupied 
state is initially characterized by the interaction energy 
t/oi and in the final state by Uq2 the transition frequency 
wiU be shifted by A[/ = C/02 - Uqi = -0.25 Uqq. Since 
L/qo is typically 1-20 kHz, this shift is sufficient to ensure 
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FIG. 2: (color online) Merging of the left target well (solid 
blue) with its right neighboring well (dashed), and elimination 
of doubly occupied wells by removal of the atom in the excited 
state conditioned on the presence of an atom in the ground 
state (a). The probability for a singly-occupied ground state 
is 1 — e, both before and after the application of the protocol. 
Wanted (unwanted) configurations are marked with a green 
crotched (red cross). This part of the protocol in itself does 
not yet increase the fidelity. 



a clean spectroscopic discrimination between the singly 
and the doubly occupied states. Following the transfer 
of one of the two atoms to another hyperfine state it is 
subsequently removed using, e.g., a resonant light pulse. 
This leaves three configurations with a single atom in the 
well. However, different vibrational levels are now occu- 
pied and they cannot be mapped to the ground state by 
a unitary process. Defining a defect as a missing popu- 
lation of the ground trap level, the error probability is 
thus still e and, so far, no improvement in the ground- 
state occupation has been achieved. The atom in the first 
excited state ((HI) in Fig. [5]), can, however, be used as 
an ancilla particle, as described in the following. 



2. Merging of the central and the left well 

In a second step, we merge the central target well with 
its left neighbor, such that the ground state of the left 
well becomes the second-excited state of the central well. 
The protocol is shown in Fig. [31 where all possible con- 
figurations and their evolutions are depicted. 

After the merging, two configurations show double oc- 
cupancy of the well, either in states \0)^ and |2)^ (I) or 
in states |1)^ and |2)^ ((HI), dashed frame). For the 
latter, we use the particle in |1)^ as ancilla to transfer 
the atom in |2)^ to |0)^ m a Sv = —2 Raman transi- 
tion, exploiting the energy shifts caused by the state |1)^ 
population (letter (a) in the figure). To see why this can 
leave the remaining configurations unchanged, we again 
need to calculate the difference in interaction energy for 
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FIG. 3: (color online) Second step of the protocol: Merging 
of the target well (solid blue) with its left neighboring well 
(dashed), which may or may not be singly occupied. The 
combination of the three distinct outcomes of the first step 
of the protocol (see the right column in Fig. [2]) with the pos- 
sible occupation or non-occupation of the left well gives rise 
to six initial configurations, two of them with an occupied 
target-well ground state (green crotched) and four without 
occupation of the target state ground state (red crosses). The 
ground state of the left well is adiabatically merged to become 
the second excited state of the target well. In (III), the atom 
in the second excited state is transferred to the ground state, 
conditioned on the presence of an atom in the first excited 
state (a). After the transfer, possible population of excited 
states is removed (b). The final transfer steps on the doubly 
occupied states result in a singly occupied ground state in 
(I)- (III). The crucial part of the protocol is thus the conver- 
sion of an initial state without occupation of the middle well 
vibrational ground state to a successful final state with unit 
occupation. 



the various configurations. For configuration (III) we ob- 
tain AU = C/oo(l - 0.875) = 0.125 C/qo, whereas the in- 
teraction energy of configuration (I) is either unchanged 
or reduced {AU < 0) if atoms are transferred to higher 
bands. 



Both two-atom configurations now populate the 
ground vibrational state with a single atom. Hereafter, 
all population in excited vibrational levels is mapped 
to the other hyperfine level {/3) via a lower sideband 
dv = —1 Raman transition and subsequently removed 
(letter (b) in the figure). This lower sideband transi- 
tion leaves the vibrational ground state population unaf- 
fected. 



Reduction of vacancies 



The protocol thus transforms the initial |1,0, 1) popu- 
lation (that occurs with probability e(l — e)^, see Tab.|l]) 
into a successful final configuration, such that three of 
the six configurations in the right column of Fig. [3] now 
have a single atom in the target well ground state. The 
probability to find a vacancy in the target well thus de- 
creases from its initial value e to 



e-e(l-e)2^2e2 



(6) 



Clearly, e' < e, and if e is small the improvement is dra- 
matic. 

The steps can be repeated iteratively, either in a par- 
allel way, such that all groups of three wells are merged 
simultaneously, or in a protocol in which only the target 
well is merged with a sequence of pairs of unpurified ad- 
jacent wells. We find the recursive update relations after 
i > 1 steps, 



- 2f2 _ f3 



2e»-i„ 



(7) 



for the cases of global parallel and sequential improve- 
ments, respectively. 



4- Unwanted transitions 

During the density-dependent transfer (step (a) in 
Fig. 131), unwanted transitions can occur, and the atom 
in the ground state, in lines (I) and (II), may be ex- 
cited. The probability for this unwanted processes de- 
pends on the pulse duration and on the detuning. For a 
driven two-level transition initially in the ground state, 
the excited-state population obeys 
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1 - cos(m)) 



(8) 



where x is the Rabi frequency, ft = \/'X^ -V is the gen- 
eralized Rabi frequency, and A is the detuning from res- 
onance. On resonance this becomes Pe = 5 (1 — cos(xi)) 
and complete transfer is obtained at t = tt/x- For the 
A = 0.125 Uoo relevant for the merging scheme (a) in 
Fig. [3] the probability of making an unwanted transfer 
to the other configurations can then be calculated using 
Eq. dl]). 

Absorption of lattice photons and collisions with atoms 
from the background gas fundamentally limit the lifetime 
r. Optimizing under the constraint r = 1 s, we obtain an 
error of £/ ~ 7 • ^-^(l • 10^*) for Uoo/h = 1(20) kHz at 
t =7(1) ms pulse duration. Thus, the lattice occupation 
error may approach, but not go below this value. Still, 
saturating this bound yields a significant improvement, 
as we will assess in the next Section. 
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C. Feasibility 

The merging algorithm can be reahzed experimentahy 
both with and without single-site resolution. The former 
is experimentally challenging but can be achieved in a 
number of different ways with a combination of existing 
technologies. First, recent advances in high resolution 
optical microscopy [1, 0] have opened the possibility to 
perform the merger either using a focussed tweezer [j| or 
with holographic [s*] or acousto-optic ^2§\ manipulation 
of the lattice potential. Without optical single-site reso- 
lution the merger can be realized by the combination of 
the flipping of the spin of an individual site (or plane) 
using a magnetic field gradient [13], spin dependent 
transport [25|, and finally double well merger in a pe- 
riod two superlattice. This approach allows merging to 
be performed on entire 2D planes in parallel. Following 
the first iteration of mergers, the neighboring sites in the 
auxiliary direction are vacant. For subsequent iterations, 
atoms have to be brought in over distances that increase 
linearly with the number of iteration steps. The analysis 
of Ref . (2^] shows that the transport fidelity over several 
lattice sites does not necessarily decrease with distance. 
Furthermore, since the reduction in state error per itera- 
tion is large, as illustrated in Eq. [71 only a few iterations 
will be needed. 

Alternatively, instead of repeated application of the fil- 
tering and merging steps discussed above, one may (adi- 
abatically) merge the ground states of n wells to the n 
lowest vibrational states of a single well in a single opera- 
tion using an appropriate n-periodicity superlattice. The 
simplest realization in a period 3 superlattice consists in 
manipulations almost identical to the ones sketched out 
in Fi g, pi Such a lattice has been implemented, e.g., in 
Ref. M. 



IV. APPLICATIONS AND ASSESSMENT 

A. Homogeneous Mott insulator 

To illustrate the accomplishments of the algorithm, we 
first apply it to a homogeneous MI state. In the zero- 
tunneling limit we calculate the initial population distri- 
bution at a given temperature using a grand canonical 
ensemble Ansatz for the number distribution: 



10" 



10" 



Q- 10 

I 



10" 



10" 




Initial distribution 



Filtering to or 1 atoms 

1 purification step 

2 purification steps 



0.05 



0.1 



0.15 
T. 



0.2 



0.25 



0.3 



FIG. 4: Defect probability (probability for vacant or multiply 
occupied sites) as a function of the dimensionless temperature 
Tjj defined in Eq. (|12p {fiu — 0.5) in the zero-tunneling MI 
regime. Initial population calculated according to Eq. (jOj. 
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(11) 
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In Fig.Ulthe solid curve shows the temperature depen- 
dence of the probability for defects (occupations other 
than unity) for a uniform MI with fiu — 0.5. With suffi- 
cient cooling this probability can be brought arbitrarily 
close to zero, but at currently achievable temperatures 
{Tu = 0.1 — 0.2 @, 0) the error probability remains ap- 
preciable. At this chemical potential, applying the fil- 
tering operation of Section |IT1 or any other mechanism 
[lOl [l3l - ll5j only slightly reduces the error probability: 
Vacancies rather than defects constitute the main source 
of error. In contrast, a single iteration of our (ideal) fil- 
tering and merging algorithm (as described in Section 
IIII[) offers more than an order of magnitude improve- 
ment as shown by the dotted curve. For example, for 
an achievable temperature of Tij = 0.1, the initial error 
probability of 10~^ is reduced to 10^^, where the largest 
contribution of the improvement stems from the merging 
algorithm. Each subsequent iteration reduces the error 
further according to Eq. 0. 



^ exp (/? [tipn - j^njn ~ 1)]) ^ 

where fj,o is the chemical potential, /3 = l/kgT is the 
inverse temperature, and 



Z = ^exp(/? 



[io-n -^n(n - 1) 



(10) 



We see that the distribution depends solely on the two 



B. Inhomogeneous Mott insulator 

In most optical lattice experiments, the atoms experi- 
ence also a weak harmonic confinement 



(13) 



In the combined potential, the atoms have an increasing 
density towards the center of the harmonic trap, with MI 
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After filtering to or 1 atoms per site (dashed line), only 
vacancies remain as an error source. The application of 
a merging step then removes the vacancies and ensures a 
uniformly filled lattice, which is extended further by the 
second merging step. A very pure domain with n = 1 is 
already formed after the first iteration and the size of this 
domain grows with subsequent iterations. Rather than 
the sharp transition of the T — MI (thick solid line), 
the pure domain is, however, surrounded by a region in 
which sites have, both, an appreciable probability to be 
occupied and to be unoccupied. In other words, this re- 
gion contains a high entropy. 



H{r) ~ — r) log P(n, r). 



(14) 



FIG. 5: Probability for unit occupation as a function of the ra- 
dius, for — 0.5. In the initial distribution (solid thin line), 
both, vacancies and multiply occupied sites are not negligi- 
ble. The filtering operation brings higher occupied sites to 
unit occupation, which increases the probability for the latter 
(dashed line). Remaining vacancies are partially eliminated 
by the application of one (dotted line) or two (dashed-dotted 
line) merging steps, which increases the probability for unit 
occupation and leads to a large close-to-ideal MI domain. 
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FIG. 6: Entropy per site, as a function of the radius, cal- 
culated according to Eq. (|14[) . While the filtering and the 
merging steps dramatically reduce the on-site entropy for low 
radii, the low average particle number per site prevents this 
for larger radii. 



domains separated by superfluid regions j|27i] . To investi- 
gate how our algorithm affects the global density distri- 
bution in this case we choose experimentally motivated 
parameters for ^^Rb Wtrap = 27r-80Hz, Uint = l.OkHz, lat- 
tice spacing d^at = 0.5/im Q. At each radius r we define 
a local chemical potential /iioc('') = /^o — Vharm('') and 
apply it to Eq. © to calculate the position-dependent 
number distribution P{n, r). 

For Tjj = 0.2 and ^jj = 0.5, we obtain the probabil- 
ity for unit occupancy plotted as the thin solid curve in 
Fig. [S] Clearly, due to the finite temperature, a signifi- 
cant part of the sites contain other than a single atom. 



which is also depicted in Fig. [51 

For applications in quantum computation this high- 
entropy shell may be of no importance as long as the 
inner region is sufficiently pure. However, for applica- 
tion in quantum simulation of spin Hamiltonians it will 
be crucial to realize a state with a high overlap with 
the full quantum many-body ground state. This can be 
achieved by using a skimmer to remove all atoms beyond 
a certain cut-off radius. Such an operation can either 
be realized with an optical tweezer focused on a single 
lattice site |28| or with magnetic single-site spectroscopic 
resolution [9|, |2J| in a cylindrically symmetric magnetic 
field configuration. We quantify the purity of the states 
of the system as the overlap (OL) with the ideal many 
body state with precisely one atom in every lattice well, 
i, of a system with radius rcut, 



0L= n ^^(1) 



(15) 



A finite value of this overlap ensures that one has a fi- 
nite component of the relevant state in, e.g., a quantum 
simulation or computation. 

In Fig. [7] we plot the infidelity, 1 — OL, versus the atom 
number in a perfect 2D MI state of size rcut ■ The upper 
four curves at = 100 illustrate how in the low filling 
case (pu = 0.5) a high-fidelity region is formed upon re- 
peated application of our algorithm. For three iterations 
a state with around 100 atoms and an overlap better than 
0.99 is realized taking into account filtering and merging 
errors (thick line), i.e. a minimal defect probability of 
ef — 1 ■ 10~^ is assumed. In the last two curves we il- 
lustrate the intrinsic power of the filtering operation at 
high filling and low enough temperature. With squares 
we show the result of the filtering on a fijj = 2 state 
with initial temperature Tjj — 0.2. As can be seen, a 
high-fidelity state of around 200 atoms is formed. Three 
iterations of our algorithm yield a nearly 10'^ atom sys- 
tem with very high fidelity: the core of an impressive 
large scale quantum computer. 

In the case of fermions, Eq. ^ contains only zero or 
one atom in the ground state but our initial filtering has 
to be extended to incorporate the possibility of zero and 
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unit population in the vibrationally excited states. Such 
filtering is simply implemented by i5i^ = — 1 Raman tran- 
sitions and removal of the final internal state atoms, 
which leaves the vibrational ground state (a) popula- 
tion intact. The merging protocol of Section IIIII relies 
on interactions, which vanish for fermions occupying the 
same hyperfine state (a) . Interactions between the differ- 
ent hyperfine states of the final state can be present and 
entropy removal by merging and occupation dependent 
filtering can be realized also in the fermionic case. 




10' 10^ 10^ 



atom, Ml 

FIG. 7: Infidelity of the states with all population beyond a 
critical radius removed shown as a function of the system size. 
Initially Tu = 0.2. The dashed line indicates the minimum 
infidelity imposed by 1-10~^ errors in the filtering and merging 
operations. 



the vibrational states occupied by the atoms. This can 
be used to filter all higher occupations to unity and thus 
remove entropy from the system. However, filtering does 
not solve the fundamental problem of vacancies, which 
remains important for low temperatures. We address 
this problem by adiabatically merging the contents of 
several lattice sites, followed by filtering, adapted to the 
vibrational states occupied after the merging. Our proce- 
dure although experimentally challenging can be realized 
with a combination of existing technology. It is based on 
global operations, and does not necessarily require single- 
site addressing, as in other proposals _IG*. As we have 
shown in Section llVAl a considerable increase in fidelity 
is achieved by the application of the merging procedure. 

The focus of our work has been on identifying effi- 
cient and fast operations, which we believe provide an 
important step towards the creation of near-zero entropy 
states of atoms in optical lattices that are useful for quan- 
tum simulations and large-scale quantum computation. 
Finally, we note that the sensitivity to atomic occupa- 
tion number using vibrational excitations may find ap- 
plications in the control of many-body quantum states 
of atoms in lattices and microtraps, e.g., by inducing en- 
ergy shifts of selected Fock components and thus tailoring 
completely new varieties of Bose-Hubbard-like Haniilto- 
nians. 



V. CONCLUSIONS 

In conclusion, we have presented an efficient way to 
realize density-dependent manipulation of optical lattice 
sites using the dependence of the interaction energy on 
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